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We perform a detailed analysis of the irreversible magnetization data of Salem-Sugui et al. and 
Babi'c et al. of underdoped and optimally doped YBa2Cu307_<5 single crystals. Near the zero field 
transition temperature we observe extended consistency with the properties of the 3D-xy universality 
class, even though the attained critical regime is limited by an inhomogeneity induced finite size 
effect. Nevertheless, as T c falls from 93.5 to 41.5 K the critical amplitude of the in-plane correlation 
length £ a bo, the anisotropy 7 = £ a f>o/£cO and the critical amplitude of the in-plane penetration depth 
Xabo increase substantially, while the critical amplitude of the c-axis correlation length £ c o does not 
change much. As a consequence, the correlation volume V^ orr increases and the critical amplitude 
of the specific heat singularity A~ decreases dramatically, while the rise of X a bo reflects the behavior 
of the zero temperature counterpart. Conversely, although £ a &o and A a (,o increase with reduced T c , 
the ratio X a bo/^~ bo , corresponding to the Ginzburg-Landau parameter K a b, decreases substantially 
and YBa2Cu307_5 crosses over from an extreme to a weak type-II superconductor. 

PACS numbers: 74.25.Bt, 74.25. Ha, 74.40.+k 



Since the discovery of cuprate superconductors, fluc- 
tuation contributions to the specific heat, magnetization, 
magnetic penetration depths etc. have been measured in 
a variety of compounds. In principle, fluctuation studies 
yield such important information as the universality class 
onto which these superconductors fall and their effective 
dimensionality^^ Given the universality class, various 
properties, e.g. the correlation volume above and be- 
low T c , arc no longer independent but related in terms 
of universal coefficients. As a consequence, the isotope 
and pressure effects on various properties arc no longer 
independent In nearly optimally doped Yh^CusOj-s, 
HgBa2Cu04_5 and La2-a;Sr :E Cu04 the occurrence of 3D- 
xy criticality is reasonably well established^^ How- 
ever, in the underdoped regime the anisotropy increases 
with reduced T c £ Although reduced dimensionality is ac- 
companied with enhanced fluctuations the 3D-xy criti- 
cal regime is expected to shrink. Nevertheless, consis- 
tency with 3D-xy scaling was observed in magnetiza- 
tion measurements of YBa2Cu307_5 down to T c ~ 61.4 
KA On the other hand, measurements of the magnetic 
penetration depth uncovered charged critical behavior 
in YBa 2 Cu 3 06.95 sample with T c ~ 61.4 K£ Note that 
7 — S = 6.95 is close to the hole concentration p ~ 1/8 
where charge fluctuations are important.- 

Recently, magnetization measurements have been per- 
formed on underdoped YBa2Cu307_a single crystals 
with T c ~ 41.5 K and T c ~ 62 K by Salem-Sugui et 
alMLH Here we perform a detailed analysis of these data. 
In contrast to previous work4^ we do not establish the 
consistency with the 3D-xy scaling plots only, but es- 
timate, given the critical exponent of the correlation 
lengths, v ~ 2/3, the critical amplitudes of the corre- 
lation length, the universal ratios, etc., of the associated 
fictitious homogeneous system as well. Indeed, the uni- 
versality class to which a given experimental system be- 
longs is not only characterized by its critical exponents 
but also by various critical point amplitude ratios and 
universal coefficients. This is achieved by invoking the 



limiting behavior of the universal scaling function, allow- 
ing to explore the growth of the in-plane and c-axis corre- 
lation lengths as T c is approached. The observed limita- 
tions of this growth are traced back to a finite size effect, 
whereupon the correlation lengths cannot grow beyond 
the extent of the homogenous domains. Clearly, such an 
analysis does not discriminate between intrinsic or extrin- 
sic inhomogeneities, but it provides lower bounds for the 
extent of the homogenous domains seen by the relevant 
fluctuations. 

The paper is organized as follows: Next we present a 
short sketch of the scaling theory and the universal prop- 
erties appropriate of anisotropic extreme type-II super- 
conductor exhibiting in the absence of an applied mag- 
netic field 3D-xy criticality. On this basis we analyze the 
irreversible magnetization data of Salem-Sugui et alMLH 
and Babic et alA We observe close to the zero field T c re- 
markable consistency with the scaling and critical proper- 
ties of a finite system belonging to the 3D-xy universality 
class. Indeed, the homogeneity of the samples turns out 
to be of finite extent, preventing the correlation lengths 
to grow beyond the extent of the homogeneous regions. 
Accordingly, the magnetization data does not provide es- 
timates of the scaling and critical properties only, but 
uncovers the spatial extent of the homogeneous regions 
as well. As T c falls from 93.5 to 41.5 K we observe that 
the critical amplitude of the in-plane correlation length 
£ a i,o, the anisotropy 7 = £, a bo/£,cQ and the critical am- 
plitude of the in-plane penetration depth \ a bo increase 
substantially, while the critical amplitude of the c-axis 
correlation length does not change much. As a conse- 
quence, the correlation volume V^ rr increases and the 
critical amplitude of the specific heat singularity A~ de- 
creases dramatically, while the rise of A a ho reflects the be- 
havior of the zero temperature counterpart^ Conversely, 
although £ a bo and X a bo increase with reduced T c , the ra- 
tio Xabo/^aboi corresponding to the Ginzburg-Landau pa- 
rameter K, a b, decreases substantially and YBa 2 Cu307_5 
crosses over from an extreme to a weak type-II supercon- 
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ductor. The rise of the anisotropy with reduced T c is con- 
sistent with a previous magnetic torque studjii^. For the 
extent of the homogenous domains we derive from the fi- 
nite size effect in £ a {, and £ c the lower bounds L a b ~ 367A, 
L c ~ 53 A and L a b ~ 254 A , L c ~ 53 A for the samples 
of Salem-Sugui et al^LLL with T c ~ 41.5 and T c ~ 62 K, 
respectively. 

To derive the scaling form of the magnetization in the 
fluctuation dominated regime we note that the scaling 
of the magnetic field is in terms of the number of flux 
quanta per correlation area. Thus, when the thermal 
fluctuations of the order parameter dominate the singular 
part of the free energy per unit volume of a homogeneous 
system scales ao 1 > 2 > l i i i 14 i li ?i ll ?i 1 7i 1 9 



where 
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This scaling form is similar to Prange's^S result for Gaus- 
sian fluctuations. More generally, the existence of the 
magnetization at T c , of the penetration depth below T c 
and of the magnetic susceptibility above T c imply the fol- 
lowing asymptotic forms of the scaling functioniSiSiiLi& 



1 dG 
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= Q c (lnz + c x ) , 
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Q is a universal constant and G (z) a universal scal- 
ing function of its argument, with G ± (z = 0) = 1. 
7 = iab/ic denotes the anisotropy, £ a t, the zero-field in- 
plane correlation length and H the magnetic field applied 
along the c-axis. Approaching T c the in-plane correlation 
length diverges as 

U = d |tr ,t = T/T c -l,± = sgn(t). (2) 

Supposing that 3D-xy fluctuations dominate the critical 
exponents are given byiS 

v ~ 0.671 ~ 2/3, a = 2v - 3 ~ -0.013, v ~ 0.671, (3) 

and there are the universal critical amplitude 
relatWiVMiW.W 



Co C ^ 



A- 



and 



A- (C 60 ) 2 Co = ^^)! = (jrj'.ff- c 0.815, 

(5) 

where A*- is the critical amplitude of the specific heat 
singularity, defined as c = {A ± /a) \t\~ a + B. Fur- 
thermore, in the 3D-xy universality class T c , and 
the critical amplitude of the in-plane penetration depth 
X a bo are not independent but related by the universal 
relation 1 ! 2 ! 14 ! 15 . 16 ! 19 , 
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From the singular part of the free energy per unit vol- 
ume given by Eq. we derive for the magnetization 
per unit volume m = M/V = ~df s /dH the scaling form 



Q ± k Bl 



TH i/2 



3/2 



F(z), F ± (z) = z 



dG 
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z dz 



with the universal coefficients^ 6 - 



o • 



-0.7, 



0.9, q = Q ± ct Q 



0.5. 



(8) 
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The scaling form J7J) with the limits (JSJ, together with 
the critical exponents (Eq. (jSJ) and the universal re- 
lations Q and |JBJ are characteristic properties of the 
3D-xy universality class. Accordingly, a homogeneous 
extreme type II superconductor falls into this universal- 
ity class when these relations are satisfied. When this is 
the case the doping dependence of the non-universal criti- 
cal properties, such as transition temperature T c , critical 
amplitudes of correlation lengths c0 , anisotropy 7, 
etc. can be determined, while the universal relations are 
independent of the doping level. 

To determine T c we consider the limit z — > 00 (x — > 0). 
Here the scaling form (J7J reduces with Eq. JSJ) to 



k B q 



H 1 / 2 



3/2 



IT, q = Q 
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(10) 



Q + c+ = Q~c^. follows from the fact that m/\/Hc adopts 
at the zero-field transition temperature T c a unique value 
where the curves m/y/H vs. T taken at different fields 
H cross and mj H X I 2 ^T C adopts the universal value 



k B q 



3/2 ' 



(11) 



Furthermore, at T c and in the limit H — > Eq. (|10|) also 
implies 



k B qj 1 

$ 3/2 #1/2 ' 



(12) 



describing the divergence of the diamagnetic susceptibil- 
ity at T c when H — ► 0. Accordingly, the location of a 
crossing point in m/yH vs. T provides an estimate for 
the 3D transition temperature and the factor of propor- 
tionality in raj (HT C ) vs. H 1 / 2 probes the anisotropy 7. 
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Given then T c and with that the reduced temperature 
t = T/T c — 1, the dominant fluctuations and their prop- 
erties can now be explored by invoking the scaling form 
0. Considering the plot M/ {TH 1 / 2 ) vs. bt for vari- 
ous fixed magnetic fields H, according to Eq. J7J) 3D-xy 
fluctuations arc verified when b scales as b cx H x l 2v with 
v ~ 2/3. As YBa2Cu307_5 near optimum doping is 
concerned, the evidence for this 3D-xy scaling behavior 
is well established^. Even though this is an essential 
step, much more detailed confirmation of 3D-xy univer- 
sality and the doping dependence of characteristic critical 
properties can be derived by invoking the limiting forms 
iJSJ of the scaling function. Considering the limit z — * 0, 
Eq. reduces below T c to 



m 
T 



~c k B 



111 



and above T r to 



m 
TH 



Cl 
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Thus, given the magnetization data of a homogenous sys- 
tem, attaining the limit z = H (^ b0 ) |i| _2l/ /<I>o << 1, 
the growth of and £ c is unlimited and estimates for 

Co. Cbo and (Cbof /Co can be deduced from 



-2/3 rn 
T 



$0 



and 



initr 4/3 + Cl 



(15) 
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given the values for the universal constants Q c and 
Q + Cq (Eq. |j5J). However, for data taken at fixed mag- 
netic field the reduction of z = (h (£+&o) 2 /*o) |£|~ 4/3 
unavoidably implies an increasing reduced temperature t 
and with that a run away from criticality. Thus, for fixed 
magnetic field the window where these limiting forms 
apply is limited to intermediate values of the reduced 
temperature. Extrinsic limitations arise in the presence 
of inhomogeneities. In this case the correlation lengths 
£,ab,c cannot grow beyond the lengths £ a &,e, set by the 
respective extent of the homogenous domains. Indeed, 
there is the Harris criterion^, which states that short- 
range correlated and uncorrelatcd disorder is irrelevant 
at the unperturbed critical point, provided that the spe- 
cific heat exponent a is negative. Since in the 3D-XY 
universality class a is negative (Eq. a rounded tran- 
sition uncovers a finite size effect^Si 2 ^, where the corre- 
lation lengths £, a b,c = Cbo cO 1*1 " cann °t grow beyond 
L a b.c, the respective extent of the homogenous domains. 
Hence, as long as £ a fc.c < -^a&.c the critical properties of 
the fictitious homogeneous system can be explored with 



the aid of Eqs. (|15fl and (|16(l . However, closer to T c the 
finite size effect sets in as £, a b,c approaches i Q 6 ;C . When 
£ c < L c and £ a & reaches L a b a finite size effect appears 
in the plot \t\~ 2 ^ m/T vs. ln|t| -4 ^ 3 around ln|t a b p | _4 ^ 3 

(£aM) \tabp\ = L a b) as the onset of deviations from 
the linear behavior. Even closer to T c where both £ a {, 
and £ c attain the respective limiting length, (m/T) tends 
according to Eq. (fLl)) to diverge as f \t\~ 2 ^ 3 where 
fo = - {Q- CQ -k B ) / ($ L C ) (ln(HL 2 ab /$ ) + Cl ). Ac- 
cordingly, sufficiently extended magnetization data do 
not provide estimates for the critical properties of the 
associated fictitious homogeneous system only, but un- 
cover the extent of the homogenous domains as well. 
As a unique size of the homogeneous domains is un- 
likely, the smallest extent will set the scale where the 
growth of the respective correlation length starts to de- 
viate from the critical behavior of a homogenous sys- 
tem. Hence, the analysis of magnetization data allows to 
probe the homogeneity of the sample as well, and pro- 
vides a lower bound for the extent of the homogenous 
domains. Furthermore, when the outlined analysis of 
magnetization data uncovers 3D-xy universality, it also 
implies that the pressure and isotope effects on the crit- 
ical properties are not independent but related in terms 
of universal relations such as 0-© and HlOfl . As an 
example, the universal form H10|) implies that the pres- 
sure or isotope exchange induced changes of magnetiza- 
tion, transition temperature and anisotropy are related 
by Am (T c ) /m (T c ) - AT C /T C + A 7 (T c ) / 7 (T C )A Last 
'but not least, it provides estimates for the pressure and 
isotope effects on the critical amplitudes, the correlation 
lengths, the transition temperature, the anisotropy and 
the extent of the homogeneous domains. 

We are now prepared to analyze the reversible magne- 
tization data of Salem-Sugui et aZ>i2*ii for underdoped 
YBa2Cu307_,5 single crystals with T c ~ 41.5 K and 
T c ~ 62 K . From magnetic torque measurements it is 
known that in the underdoped regime the chemical sub- 
stitution tuned reduction of T c is accompanied by an in- 
crease of the anisotropy. From Fig. ^ showing 7 vs. T c 
it is seen that 7 increases from 6 around T c ~ 91 K to 29 
near T c ~ 40 K. 

It is well documented that in highly anisotropic 
cuprate superconductors the magnetization curves 
M(T) taken at varying fields cross below T c around 
(T cross ,M cross ), nearly independent of the magnitude 
of the applied field24. Although this crossing phe- 
nomenon is not associated with criticality, it uncovers 
that the system behaves below r cross nearly quasi two 
dimensionaiiaSi 2 ^ In Fig. [5] we depicted the data 
of Salem-Sugui et aZiISiil for underdoped YBa2Cu307_^ 
single crystals with T c ~ 41.5 K and T c ~ 62 K 
in terms of M vs. T for varying H applied along 
the c-axis. Apparently there are nearly field indepen- 
dent crossing points around T cross ~ 40.3 K and 60.5 
K. This phenomenon is however absent in nearly opti- 
mally doped and with that considerably less anisotropic 
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FIG. 1: 7 vs. T c for YBa2Cu307-,5 derived from Janossy 
el alw— To estimate 7 for the T c 's considered here we use 
7 = 1150/r c indicated by the solid line. 



YBa2Cu307_,5 samples^ Since 7 increases with reduced 
T c (see Fig. the occurrence of a nearly field inde- 
pendent crossing point uncovers then a 3D to quasi 2D 
crossover below T cross . 

Nevertheless, above T cross 3D-xy fluctuations are ex- 
pected to set in and to dominate. In this case Eq. ljTU|) 
implies the occurrence of a crossing point in M/(TH 1 / 2 ) 
vs. T at T c . The plots shown in Fig. uncover these 
crossing points and provide for the respective transition 
temperatures the estimates T c ~ 41.5 K and T c ~ 62 K. 

To identify the dominant fluctuations wc rcplottcd the 
data depicted in Fig. Ok in terms of M / (TH 1 / 2 ) vs. bt 
with b adjusted to achieve a collapse on the H = 0.1 T 
curve close to t = 0. According to Eq. JIJ) the domi- 
nance of 3D-xy fluctuations is established when b scales 
as b oc l/H 1 / 2 " oc I/H 3 / 4 . In Fig. Hd, depicting the field 
dependence of b, we observe remarkable consistency with 
the characteristic 3D-xy behavior. However, the quality 
of the data collapse is seen to deteriorate with increasing 
field. 

According to Fig. [5] showing the corresponding plots 
for the sample with T c ~ 62 K, we observe again con- 
sistency with the characteristic 3D-xy critical behavior, 
b ~ H~ 3 / 4 . Even though the data collapse is seen to 
deteriorate with increasing field as well, these plots un- 
cover the dominance of 3D-xy fluctuations around the 
estimated transition temperature T c . 

Having established the consistency with 3D-xy univer- 
sality in terms of scaling plots based on the universal 
scaling form J7J) with v ~ 2/3 and the estimates for T c 
we are no prepared to perform a more detailed analysis 
providing additional checks, allowing to explore the inho- 
mogeneity induced finite size effects and to estimate the 
essential critical properties, including the critical ampli- 
tudes of the correlation lengths and the anisotropy. The 
basic starting point is either the limiting behavior below 
or above T c , given by Eas. (|15|l and i|16|) . respectively 
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FIG. 2: Reversible magnetization M vs. T for various fixed 
magnetic fields H applied along the c-axis of underdoped 
YBa2Cu307_<s single crystals with T c ~ 41.5 K (a) and 
T c ~ 62 K (b) derived from the data of Salem-Sugui et alJ&l± 
There are nearly field independent crossing points around 
T cross ~ 40.3 K and 60.5 K. 



TABLE I: T c , weight of the samples and relationship between 
M in emu and m = Mp/weight with p ~ 6.3 g/cm 3 . 



T c (K) 


weight (mg) 


m (emu cm 3 ) 


41.5 


1 


6300 M 


62 


1.2 


5250 M 



To invoke these limiting forms it is necessary to convert 
the magnetization data given in emu to m in emu cm~ 3 
according to Table IJ 

Next we invoke the limiting form Eq. I|15|l to estimate 
the magnitude of the critical amplitudes below T c and to 
explore the homogeneity of the samples. For this purpose 
we depicted in Fig. \t\~ 2/3 m/T vs. ln|tf 4/3 . In both 
samples we observe in a limited interval consistency with 
the asymptotic behavior, indicated by the straight lines. 
Using Eqs. © and l|15|) we derive from the slope of these 
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FIG. 3: M/H 1/2 vs. T at various fixed magnetic fields H 
applied along the c-axis of underdoped YBa2Cu307_<5 single 
crystals with T c ~ 41.5 K (a) and T c ~ 62 K (b) derived 
from the data of Salem-Sugui et aZ^fi*— The crossing points 
provide an estimate for the respective transition temperature, 
T c ~ 41.5 K (a) and T c ~ 62 K (b) 



lines the estimates 

Co =i 1.46 A (T c ~ 41.5 K), Co - 1-33 A (T c ~ 62 K), 

(17) 

for the critical amplitude of the c-axis correlation length. 
From the straight lines in Fig. it also follows that 

10 3 (r ) 2 

In = -5.09 - a (T c ~ 41.5*0, 



In 2-5 10 ' ^ g6o) = -5.39 - c x (T c ~ 62), (18) 

yielding for the critical amplitudes of the in-plane corre- 
lation lengths the ratio 



Cam (r c ^41.5K) ^ i5i 
C b0 * 62 K) 



(19) 



In addition, this plot reveals that the attainable criti- 
cal regime, indicated by the straight lines, ends around 
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H CO 

FIG. 4: a) Mj (tH 1/2 ^J vs. bt with b adjusted to achieve a 

collapse on the H — 0.1 T curve for the data shown in Fig. 
0i. b) b vs. H. The solid line is b oc H 3 / 4 characteristic for 
3D-xy thermal fluctuations. 



-4/3 



4.12 and In |U ' 4/3 



the marked values In \t a bp\ — auu in |iai>p| 
4.46. Here the in-plane correlation length £ ab reaches the 
limiting length L ab , so that Co \t abp \~ 2 ^ 3 = L ab . From 
In |t Q 6 p r 4/3 ^ 4.12 and In \t abp \~ i/3 ~ 4.46 and Eq. O 
we obtain for the ratio between the limiting lengths L a h 
the estimate 



Lab (T c ~ 41.5 K) 
i ah (T c ~ 62 K) 



1.36. 



(20) 



Nevertheless, given the estimate for the universal coeffi- 
cient ci , we can extract from Eq. (|18fl the critical ampli- 
tude of the in-plane correlation length of the respective 
fictitious homogeneous system. Using c\ = 1.76, which 
will be derived later on, we obtain 

Co - 46.82 A: T c ~ 41.5 K, Co - 27.32A: T c ~ 62 K 

(21) 

and together with the values for £ c0 (Eq. (JTJJ) for the 
anisotropy 



7 ~ 32.07: T c ~ 41.5 K, 7 ~ 20.5: T c ~ 62 K, 



(22) 
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FIG. 5: a) M/ (TH 1/2 ^j vs. bt with b adjusted to achieve a 

collapse on the H = 0.25 T curve for the data shown in Fig. 
|3>. b) b vs. H . The solid line is b oc H~ 3 ^ characteristic for 
3D-xy thermal fluctuations. 



in reasonable agreement with the estimates derived from 
the magnetic torque measurements, namely 7 ~ 28 and 
19 (see Fig. [TJ. 

To analyze the behavior closer to T c we consider the 
plot — \t\~ 2 ^ 3 m/T vs. —t shown in Fig. [7| Above 
t abp ~ 0.046 (Fig. EH) and t abp ~ 0.034 (Fig. Qd) there 

is an interval exhibiting the In \ t\ 4 ^ 3 behavior of a ho- 
mogeneous system, indicated by the solid curve and con- 
sistent with the \t\ 2 ^ 3 m/T vs. In | i | 4 ^ 3 plots shown in 
Fig. Indeed, a minimum occurs below t abp followed 
by an increase. Although the data is sparse, it indicates 
a \t\ 2 ^ 3 divergence, associated with a finite size effect in 
both, £ Q t, and £ c , whereby these correlation lengths can- 
not grow beyond L ab and L c , respectively. According to 
Eq. l(TH)l the amplitude of the divergence is given by 



, 1—2/3 TTl 
\t\ 2/3 y 



Q~ c o k B f ln (HLl b 
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FIG. 6: \t\ 



-2/3 < 



/T vs. In \t\ 4/3 for t < 0. a) Sample with 



-4/3 



The arrow marks In \t„ 



T c ~ 41.5 K at .ff = 0.1T; the solid line is \t\~ 2/3 m/T = 
-0.163+0.032 lnltl 



-4/3 



4.12 



(t abp ~ -0.045). b) Sample with T c ~ 62 K at H = 0.25T; 
the solid line is \t\~ 2/3 m/T = -0.189 + 0.035 In \t\~ 4/a . The 
arrow marks In \t abp \~ A ^ 3 — 4.46 (t abp ~ —0.035). 



and allows to determine i c , given the amplitude /o, L a t 
and ci. From ln|U P |" 4/3 4.12 (Fig. Et), £~ 60 - 46 - 82 
A (Eq. d)) and ]n\t abp \~ A/3 ~ 4.46, Qo - 27 - 32 A 
we obtain with £ afc (i abp ) = £~ 6 liat.^, | ~ 2/3 = L ab for the 
limiting length in the afc-planc the estimate, 

L ab ~ 367 A: T r ~ 41.5 K, L ab ~ 254 A: T c ~ 62 K, 

(24) 

in reasonable agreement with the estimate for their ratio 
given by Eq. Together with Q~c ~ -0.7 (Eq. 

©), Ci = 1.76 and the rather crude estimates for / (see 
Fig. O) Eq. 123) yields 

L c ~ 57 A: T c ~ 41.5 K, L c ~ 53 A: T c ~ 62 K, (25) 

in comparison with the previous estimates, L ab ~ 392 A 
and L c ~ 52 A for YBa 2 Cu 3 6 .7 with T c ~ 59.7 K, de- 
rived from the inhomogencity induced finite size effect in 
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the temperature dependence of the in-plane penetration 
depth close to criticalityiSi 
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FIG. 7: -|i|~ 2/3 m/T vs. -t for t < 0. a) Sample with 
T c ~ 41.5 K at H = 0.1T; the solid curve is - |t|~ 2/3 m/T = 
0.15 - 0.03 In |tr 4/3 and the dashed one - \t\~ 2/3 m/T = 8 
10~ 4 |f|~ 2/3 . The arrow marks t abp ~ -0.047 (In \t abp \~ i,z ~ 
4.12). b) Sample with T c ~ 62 K at H = 0.25 T; the solid 
line is - \t\~ 2/3 m/T = 0.19 - 0.035 In |t|" 4/3 and th e dashed 
one - |t|~ 2/3 m/T = 0.0008 |i|~ 2/3 . The arrow marks t abp ~ 
-0.042 (ln|t abp r 4/3 -4.46). 



Further estimates of the critical amplitudes of the cor- 
relation and limiting lengths can be obtained by invoking 
Eq. (O- According to this we plotted |t| 2/3 m/ (TH) 
vs. t for the sample with T c ~ 41.5 K and H = 0.2 T in 
Fig. IHt- As t tends to zero the data approaches a mini- 
mum around t = t cp ~ 0.026, indicated by the horizontal 
line and marked by the arrow. This minimum provides 
an estimate for the critical behavior given by Eq. (|14f> . 
while the upturn below uncovers again an inhomogene- 
ity induced limiting length along the c-axis and in the 
afe-plane. Indeed, when the growth of £ a b,c are limited 
by L abiC , \t\ 2/3 m/ (TH) tends to zero because the ratio 
£ab/£c approaches the ratio L 2 ab /L c . A glance to Fig. [SJ) 
reveals that the minimum exhibits a linear field depen- 



dence. Accordingly, the limiting behavior dG/dz = c^z 
(Eq. (JSJ) is not attained and the lowest magnetic field 
dependent correction, compatible with the linear depen- 
dence, dG/dz = Cq"z(1 + gz), must be taken into ac- 
count. The linear extrapolation yields 



2.08 10" 8 (emu cm _3 K _1 G _x ), 



(26) 

and with the univcr sal coefficient Q+c^ ~ 0.9 (Eqs.©) 

(C h o) 2 /e c + o = 7C6o^718A. (27) 

Invoking then our estimate for the anisotropy, 7 ~ 32.07 
(Eq. (J23)) we obtain 

Cm * 22 - 4 A > Co = Cwh * 0.7 A, (28) 
and with £^ ~ 1.46 A (Eq. I|19[l) for the universal ratio 



Co 
Co 



SafcO 



2.1, 



(29) 



compared to the theoretical prediction £,co/£cO = 
%abo/€abo — 2 - 21 - (^Q- ©)■ Unfortunately, an equivalent 
analysis of the data for the less underdoped sample with 
T c ~ 62 K is not opportune, because the measurements 
do not extend to comparably low fields. 

In Table[n]wc summarized our estimates for the critical 
properties of underdoped YBa2Cii307_5 single crystals 
derived from the magnetization data of Salem-Sugui et 
alrliLl^. For comparison we included corresponding val- 
ues for nearly optimally doped YBa2Cu3C>7_#i£. While 
the critical amplitude of the c-axis correlation length 
exhibits a rather weak doping dependence the ab- 
plane counterpart, ^ h0 , and the anisotropy increase dras- 
tically with reduced T c . As a consequence, the cor- 
relation volume V (Mrr = (£^0) £<T increases and the 
critical amplitude of the specific heat singularity A~ 
(Eq. JSJ) decreases dramatically. This behavior ren- 
ders it difficult to extract in the underdoped regime crit- 
ical behavior from specific heat data. Because £,~ does 
not change much with reduced T c , the critical ampli- 
tude of the in-plane correlation length A a bOj resulting 
from Eq. rises considerably, reflecting the behav- 

ior of its zero temperature counterpart A a j, (0), where 
A a6 (0) ~ 1300 A and A a6 (0) ~ 2250 A at T c = 91.3 K 
and 60.5 K, respectivcljii^. Conversely, although both, 
the critical amplitudes of the in-plane correlation length 
£abo ano - penetration depth A a bo increase with reduced 
T c , the ratio \ a bo/£abO' corresponding to the Ginzburg- 
Landau parameter K a b, decreases substantially, where- 
upon YBa2Cu307_^ crosses over from an extreme to a 
weak type-II superconductor. Noting again that the size 
of the homogeneous domains is not necessarily unique, 
the onset of the finite size effect probes their respective 
smallest extent. Accordingly our estimates for L a b and 
L c are lower bounds for the extent of the homogenous 
domains. 
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TABLE II: Collection of the estimates derived from the magnetization data of Salem-Sugui et alw—— for the underdoped 
YBa2Cu307-i single crystals. For comparison we included the estimates for nearly optimally doped YBa2Cu307_<5^- The 
bracketed values for and fj" 60 are obtained with Eq. J1J. V c ~ OTr = (£~ M ) denotes the correlation length volume below 
T c . To obtain the critical amplitude of the magnetic in-plane penetration depths A a f,o we used the universal relation @. 



T c 


Co 


£abO 


7 — Cbo/£cO 


c+ 


(Cbo) /5c0 


c+ 

SabO 


v- 


Lab 


L c 


10 3 A a bo 


Aabo/CabO 


10 6 A^bo T c 
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A 


A 
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A 3 


A 


A 


A 




A 2 K 


91.7 


1.3 
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0.941 
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81.2 


62 


1.33 


27.32 


20.5 


(0.60) 




(12.36) 
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254 


53 


1.157 


42 
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41.5 
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46.82 
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FIG. 8: a) \t\ 2/3 m/ (TH) vs. t for the sample with T c ~ 41.5 
K and H = 0.2 T. The horizontal line marks the minimum, 
b) \t\ 2/3 m/ (TH)\ min vs. H for the sample with T c ~ 62 



K and H — 0.25 T. The line is \t 
10" 8 + 2.3 10~ 8 H with H in T. 



2/3 , 



i/(TH)\ 



= -2.08 



In this context it is important to recognize that the 
listed anisotropics refer to the homogeneous counter- 
parts. However, in the actual inhomogeneous samples 
these values only apply in an intermediate temperature 
regime where the growth of the correlation length is not 
yet limited by the finite size effect. Indeed, consider- 
ing the sample with T c = 41.5 K, below T c £ a b lev- 
els off around \t a bp\ — 0.045, while £ c saturates around 



\t cp \ ~ 0.0041. Accordingly, below \t a b P \ the anisotropy 
7 = iab/ic decreases from 7 = Q> /£tf) - 32 to 7 = 
Lab/L c ~ 6.4 at T c . This behavior implies that the 2D- 
limit in the underdoped regime is hardly accessible be- 
cause £ a b cannot grow beyond L a b and as a result 7 does 
not diverge for fixed £ c . 

To check the consistency of our analysis further, we 
invoke Eq. J5J to calculate from the magnetization data 
the derivative of the universal scaling function in terms 
of 



.dG 



m $ 



Q dz ~ Tk B ^° 



c ± 1.1—2/3 



= (H(^bo) 2 /^ 



-4/3 



. (30) 

and the respective estimates for the critical amplitudes 
of the correlation lengths (Table [H}. In Fig. EJi, showing 
—Q~ dG/dz vs. ln(z) below T c , we observe from ln(z) ~ 
—2.75 down to —4.25 consistency with the leading z — > 
behavior Q~dG/dz — Q~Cq (In (z) + c±), where Q~Cq — 
—0.7 (Eq. ©) and C\ was chosen as 



Cl 



1.76. 



(31) 



This estimate fixes the so far unknown universal coef- 
ficient ci.The systematic deviations, setting in around 
ln(z) ~ —2.75 (z ~ 0.065), uncover the onset of the inho- 
mogencity induced finite size effect, limiting the growth 
of £ a b. Indeed, this value corresponds to z = HL 2 b /& 
with H = 0.1T and L ab = 367 A (TableHU). On the other 
hand, the upturn setting in around ln(z) ~ —4.2 (z ~ 
0.0136) signals the escape from the scaling regime where 
the ln(z) behavior applies. Indeed, in the present case H 

is fixed and the reduction of z = (^H (£,^ b0 ) 2 /$oj |*| 4 ^ 3 
unavoidably implies an increasing reduced temperature 
t and with that a run away from criticality. Thus, 
for fixed magnetic field the window where the univer- 
sal scaling function can be observed is limited from be- 
low by the run away from criticality and from above 
by the finite size effect in In Fig. |5Jd we de- 

picted Q+d 2 G/dz 2 = -d(m^- Q \t\~ m $ / (*bT)) /dz 

vs. z = (#fctfco) 2 /$o) |*P 4/3 - According to Eq. © it 
approaches in a homogenous system in the limit z — > 
the universal value Q + d 2 G/dz 2 = Q + Cq ~ 0.9. Even 
though the available data are rather sparse we observe 
in the interval 0.01 < z < 0.04 consistency with this 
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limiting behavior, indicated by the horizontal line. Un- 
fortunately, the available data do not allow to locate the 
onset of the finite size effect in £ a (,, seen below T c around 
ln(z) = hx(HLl b /^ ) ~ In (0.065) ~ -2.73(Fig. Hi). 
The upturn setting in around z = 0.01 signals the run- 
away from criticality as in FigEti. 




o 




0.08 



FIG. 9: a) Q'dG/dz 



(m$ /Tfcs)Col* 



-2/3 



for the sample with T c = 41.5 K at H = 0.1 T, Q, = 46.82 
A and £~o = 1.46 A (O) and the sample with T c = 62 
K at H = 0.25 T, = 27.32 A and Co = 1.33 A 

(□); the solid line is Q~dG/dz = Q~Cq (In (z) + ci) with 

* = (h (Cm>) 2 /*o) |t|~ 4/3 , Q-Co = -0.7 (Eq. ®), and 
ci = 1.76 (Eq. 1311 ). The arrows mark the onset of 
the finite size effect in C at ln(z) = HL^ b /^o, namely 
Ih(jz) = -2.73 (T c = 41.5 K) and ln(jz) = -2.55 (T c = 62 

K). b) Q+d 2 G/dz 2 = -d(m($> /Tk B )Co |t|" 2/3 ) Jdz vs. 

* = {H {CoY /*o) It!" 473 with T c = 41.5 K for £+o = 22.4 
A and Co = 0.7 A; the solid line is Q + d 2 G/dz 2 = Q + cJ with 
Q+eJ ~ 0.9 (Eq. ©). The arrow marks z = HL 2 ab /$ ^ 
0.065, the onset of the finite size effect in C- 



the data of the Yl^O^O^-a single crystal with T c ~ 
93.5 K. In Fig. EH we show M vs. ln(-ff). In a limited 
interval we observe linear behavior so that the scaling 
form (|13[) rewritten in the form 



M 



VQ~ 



-0 K B 



= d 



eln(ff), 



In (H) + In 



<2 
Sab 

$0 



Cl 



(32) 



applies. The solid lines arc this scaling form with the 
parameters d and e listed in Tabic ITTT1 
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FIG. 10: M vs. ln(_ff) for various temperatures derived from 
the data of Babi'c et aii- for a YBa2Cu307_5 single crystal 
with T c — 93.5 K and the magnetic field applied along the c- 
axis. The solid lines are Eq. 13211 with the parameters listed 
in Table ITTT1 

Given the volume of the sample, the universal ampli- 
tudes Q~Cq ~ -0.9 (Eq. ©) and ci (Eq. (gU) the listed 
correlation lengths are then readily calculated. Together 




with C 



bO,cO 



W 



-2/3 



we obtain for the critical am- 



plitudes the estimates 



Cbo - 7 A, Co * 1-4 A, (33) 

in comparison with £~ b0 ~ 10.4 A and Co — 1-3 A for the 
sample with T c = 91.7 K (Table EJ. 



TABLE III: Parameters entering Eq. 
straight lines in Fig. 1111 and with Q~c 
ci = 1.76 (Eq. and V = 8.2 10~ 4 c 

for the correlation lengths C and C 



132H . yielding the 
= -0.7 (Eq. @), 
£ the estimates 



T(K) 


d (emu) 


e (emu) 


U(A) 


Cc(A) 


91 


-2.36 10~ 3 


2.13 10~ 4 


74.1 


16.3 


90 


-3.10 10" 3 


2.70 10" 4 


59.5 


12.8 


89 


-3.93 10~ 3 


3.32 10~ 4 


50.7 


10.3 


88 


-4.72 10" 3 


3.93 10" 4 


46.5 


8.8 



Next we turn to magnetization data of Babic et ai£, 
taken at fixed temperatures below T c as a function of the 
magnetic field applied along the c-axis. Here we analyze 



To check the consistency with the scaling form <|32[) 
further, we calculated m$o£ c / (ksT) = —Q~dG/dz vs. 
z for T = 91 K and 90 K as shown in Fig. El 
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with the parameters listed in Table II I II The compar- 
ison with the leading z — > behavior, —Q~ dG/dz = 
—Q~Cq (ln(z + ci), reveals that this regime is attained, 
but limited by irreversibility in the limit z — > and the 
crossover to the large z-limit , —Q~dG/dz = —qz 1 / 2 with 
q ~ 0.5 (Eq. ©). This limitation is also apparent in Fig. 
ITUl on the low and high field side. Indeed, because the 
data do not extend close to T c , the correlation lengths are 
comparatively small and their growths is not yet limited 
by the extent of the homogenous domains. 




U I ■ '■ 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



z 

FIG. 11: m$o£c/( fe sT) = -Q~ dG/dz vs. z for T = 91 K 
and 90 K. The solid line is —Q~ dG/dz — —Q~Cq (\n(z + ci) 
with Q- C g = -0.7 (Eq. ©), ci = 1.76 (Eq. (ISTIh 

Although the analyzed magnetization data is afflicted 
with uncertainties arising from the subtraction of the 
normal state paramagnetism and the Curie term due 
to paramagnetic impurities or defects, we observed re- 
markable consistency with 3D-xy critical behavior for 
both, nearly optimally doped and underdoped samples. 
In contrast to previous work^ we did not establish the 
consistency with the 3D-xy scaling plots only, but es- 
timated, given the critical exponent of the correlation 
lengths, v ~ 2/3, the critical amplitudes of the cor- 
relation length, the universal ratios, etc., of the asso- 
ciated fictitious homogeneous system as well. Indeed, 
the universality class to which a given experimental sys- 



tem belongs is not only characterized by its critical ex- 
ponents but also by various critical point amplitude ra- 
tios and universal coefficients. This has been achieved 
by invoking the limiting behavior of the universal scal- 
ing function dG/dz, allowing via Eqs. I|13|) - (|16fl to ex- 
plore the growth of the in-plane and c-axis correlation 
lengths as T c is approached from below or above. We 
have seen that this growth is limited due to inhomo- 
geneities and that this limitation appears to be equiv- 
alent to a finite size effect, whereupon the correlation 
lengths cannot grow beyond the extent of the homoge- 
nous domains. Clearly, such an analysis docs not discrim- 
inate between intrinsic or extrinsic inhomogeneities, but 
it provides lower bounds for the extent of the homoge- 
nous domains seen by the relevant fluctuations. Even 
though since the discovery of superconductivity in the 
cupratcs by Bednorz and Mullei^ a tremendous amount 
of work has been devoted to their characterization, the 
issue of inhomogeneities and their characterization is still 
a controversial issue. There is neutron spectroscopic evi- 
dence for nanoscale cluster formation and percolative su- 
perconductivity in various cuprates22^£. Nanoscale spa- 
tial variations in the electronic characteristics have been 
observed in underdoped I^S^CaC^Os+s with scan- 
ning tunneling microscopji2ii22i2iiM. They reveal a spa- 
tial segregation of the electronic structure into 3-nm- 
diameter superconducting domains in an electronically 
distinct background. Furthermore, the investigations of 
Gauzzi et alm^ on YBa2Cu306.9 films with reduced long- 
range structural order clearly reveals that the size of the 
homogeneous domains strongly depends on the growth 
conditions. In any case we have shown that the analy- 
sis of reversible magnetization data taken near criticality 
does not uncover the critical properties of the associated 
fictitious homogenous and infinite system only, but pro- 
vides lower bounds for the extent of the homogeneous 
domains as well. Last but not least, having established 
the consistency with 3D-xy universality there are univer- 
sal relations such as and (|10f) . They imply that 
the effect of pressure and isotope exchange on the respec- 
tive properties are not independent. 

The author is grateful to S. Salem-Sugui Jr. and J. R. 
Cooper for providing the magnetization data. 
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